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' The  results  are  presented  of  a combined  experimental  and  theoretical 
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20.  attributed  to  an  effect  of  excitation  in  the  NO^  beams  of  the  merged  beams 

and  it  is  concluded  that  because  of  these  effects  the  flowing  afterglow  is 

presently  better  for  obtaining  rate  coefficients  for  complex  molecular  ions. 

Our  theoretical  work  has  focused  on  two  neutralization  mechanisms: 
electron  transfer  and  tidal  distortion.  A new  formulation  of  the  collision 
dynamics,  Semlclasslcal  Perturbation  Scattering  Theory,  was  developed  and 
used  to  estimate  the  importance  of  the  tidal  distortion  mechanism.  It  is 
concluded  that  the  electron  transfer  mechanism  dominates  in  reactions  of 
simple  ions  and  small  clusters,  but  that  the  tidal  distortion  becomes  more 
important  for  multiply  hydrated  reactants.  An  approximate  scaling  formula 
for  the  rate  constants  is  obtained  that  fits  a wide  range  of  experimental 
data  involving  simple  ions  and  small  clusters  with  an  accuracy  of  about  ± 307.. 
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PREFACE 


This  Final  Report  covers  work  done  under  contract  F19628-75-C-0050, 
which  has  been  a combined  theoretical  and  experimental  study  of  mutual 
neutralization  reactions  between  positive  and  negative  ions.  These 
reactions  are  responsible  for  the  final  removal  of  ions,  terminating  a 
long  series  of  ion-neutral  reactions  that  follow  the  initial  ioniza- 
tion and  subsequent  negative  ion  formation  in  the  D-  and  sub-D  regions 
of  the  earth's  ionosphere.  Thus  the^r  must  be  included  in  modeling  the 
ion  chemistry  of  the  ionosphere,  an<f  they  play  a particularly  sensitive 
role  in  those  models  devoted  to  VLF  and  EIF  communications  problems. 
Unfortunately,  those  reaction  rates  pertinent  to  the  Vlf  and  ELF 
problems  are  not  yet  determined  with  the  desired  accuracy  and 
reliability;  however,  we  have  determined  that  the  existing  experimental 
facilities  at  SRI  are  not  capable  of  substantial  improvement  toward 
that  end. 

Because  this  report  covers  the  completion  of  the  SRI  work 'using 
the  merged-beam  technique,  and  no  further  experimental  work  is  contem- 
plated, we  take  this  opportunity  to  review  and  sunmarize  all  experi- 
mental and  theoretical  work  done  to  date,  including  experiments  at 
other  laboratories.  This  report  thus  serves  as  an  up-to-date  compilation 
and  review  of  results  on  two-body  reaction  rates.  Three* body  rates  are 
briefly  discussed  because  of  their  Importance  at  low  altitudes. 
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INTRODUCTION 


In  this  report,  we  will  discuss  recent  advances  in  the  measurement 

and  calculation  of  reaction  rate  coefficients  for  the  positive-ion— 

negative-ion  mutual  neutralization  process 

+ - 

A + B - A + B , (1) 

which  has  important  effects  on  atmospheric  ionization.  In  reaction  (1), 
A and  B may  be  molecular  ions,  and  the  neutral  products  may  be 
excited,  dissociated,  or  rearranged;  very  little  is  known  about  the 
final  products.  The  corresponding  three-body  mutual  neutralization 
reaction, 

+ - 

A + B + M - (A+B)  + M (2) 

is  important  at  gas  pressures  greater  than  a few  torr  (at  altitudes  of 
less  than  40  km). 

Eq.  (1)  is  of  current  interest  in  the  chemical  modeling  of  the 
ionosphere,  particularly  in  the  D-region,  because  it  is  the  final  step 
that  removes  ionization.  The  three-body  reaction  (2)  is  important  in 
the  stratosphere.  Three-body  mutual  neutralization  has  long  been 
studied  in  regard  to  gaseous  electronics  (see  Thomson  and  Rutherford, 

1896),  has  been  reviewed  by  Mahan  (1973),  and  has  seen  a recent  revival 
of  interest  because  of  its  role  in  the  kinetics  of  rare-gas/halogen 
laser  media  (see  for  example,  Flannery  and  Yang,  1978;  and  Wadehra  and 
Bards  ley,  1978). 

Mbtual  neutralisation  reactions  are  important  in  ionized  gases 
because  of  their  large  cross  sections,  attributable  to  the  coulomb 
interaction.  Binary  mutual  neutralization  cross  sections  at  an 
average  relative  energy  E corresponding  to  300  K are  typically 
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10  cm  and  scale  as  E at  low  energies.  Gas  phase  reaction  rate 

coefficients  for  simple  atmospheric  ions  in  the  temperature  range 

-7  3 -It 

T “ 200-300  K thus  tend  to  be  about  10  cm  /sec  and  scale  as  T for 
temperatures  T < 1,000  K. 

Despite  the  large  cross  sections,  mutual  neutralization  reactions 
are  difficult  to  study  in  the  laboratory  because  of  the  difficulty  in 
creating  sufficiently  high  densities  of  identified  and  selected  species 
of  ions.  Among  experiments  on  ternary  mutual  neutralization,  identifi- 
cation of  the  ion  species  involved  is  in  question.  Binary  mutual 
neutralization  experiments  are  in  better  shape,  but  even  so,  almost  no 


work  has  been  done  on  identifying  the  neutral  products  and  excited 
states  resulting  from  the  reactions. 

Ion-ion  mutual  neutralization  forms  an  interesting  theoretical 
problem  because  it  belongs  to  a large  class  of  electron- transfer 
processes  that  can  be  described  in  terms  of  potential  curve  crossing 
dynamics;  further,  it  represents  a case  in  which  the  shapes  of  the 
potential  curves  are  well  known.  However,  the  calculation  of  mutual 
neutralization  cross  sections  is  complicated  by  the  large  number  of 
potential  curve  crossings  involved  in  many  cases,  especially  with 
molecular  species  where  rotational  or  vibrational  excitation  may  play 
a role  in  the  collision. 


A "typical"  atomic  ion-ion  mutual  neutralization  cross 
section  as  a function  of  relative  energy 


(1)  the  large  cross  section  values,  which  are  typically  > 10,000  A 

at  an  er jrgy  corresponding  to  300  K;  (2)  an  E 1 energy  dependence  at 

very  low  energies  (E  ^ 0.1  eV);  (3)  an  E behavior  at  very  high 

_2 

energies,  which  is  really  a v velocity  dependence  (independent  of 
mass),  and  which  may  be  different  if  dissociation  occurs--see  Moseley, 
Olson,  and  Peterson,  (1975),  pp.  33-34;  and  (4)  a relatively  flat 


portion  of  the  curve,  which  may  even  increase  with  E at  intermediate 
energies.  The  last  feature  is  of  greatest  complexity  theoretically, 
because  in  that  energy  region,  the  different  potential  curve  crossings 
exhibit  the  greatest  influence.  The  cross  sections  here  are  also 
really  velocity  dependent,  thus  the  extent  in  energy  and  general 
behavior  of  the  cross  section  in  the  region  vary  considerably  among 
different  reactants.  The  cross  sections  may  have  local  maxima  and 
minima  as  various  final  states  or  reaction  channels  become  important. 

In  the  H + H system,  some  closely-spaced  undulations  in  the  cross 
sections  were  discovered  by  Moseley,  Aberth,  and  Peterson  (1970),  and 
were  interpreted  as  the  result  of  interference  between  different  channels 
For  molecules,  the  energy  dependence  characteristics  (features  2 and  4) 
tend  to  be  less  distinct;  the  intermediate  part  of  the  curve  is  both 
more  extensive  and  featureless  because  of  the  increased  number  of 
final  states. 


II  EXPERIMENTAL  TECHNIQUES  AND  RESULTS 


Several  careful  experiments  on  ternary  ion-ion  mutual  neutrali- 
zation were  conducted  in  the  1930s  in  x-ray  excited  gases,  following 
development  of  the  technique  over  many  years.  The  experiments  showed 
that  the  effective  two-body  rate  coefficient  increased  as  the  gas 

pressure  was  increased — owing  to  three-body  effects — but  saturated 
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for  gas  pressure  of  about  1 atm,  at  a value  of  2-3  x 10  cm  /sec. 
Figure  2,  for  measurements  in  air,  shows  results  of  Sayers  (1938)  as 
open  circles,  and  of  Machler  (1936)  as  solid  dots,  compared  to  the 
theoretical  model  of  Natanson  (1959).  It  is  interesting  that  despite 
the  uncertainty  in  the  identity  of  the  ionic  species  involved. 


Figure  2.  Typical  results  for 
three-body  ion-ion  mutual  neutral! 
zation  in  air,  expressed  in  terms 
of  an  effective  binary  rate  coef- 
ficient. Figure  from  Natanson 
(1959). 


saturation  values  of  2-3  x 10  cm  /sec  at  300  K are  consistently 
reported,  and  that  range  of  values  may  hold  even  for  atomic  ions 
(Wadehra  and  Bardsley,  1978). 

More  recent  gas  phase  measurements  of  Mahan,  and  of  Hirsh  and 
Eisner,  have  been  discussed  in  the  reviews  of  Mahan  (1971)  and  of 


Moseley,  Olson,  and  Peterson  (1975),  respectively.  The  interpretation 
of  the  data  is  complex. 

B.  Merged  Beams-Technique 

In  the  merged-beams  technique,  two  mass-selected  ion  beams  moving 
in  the  same  direction  are  merged  and  superimposed  over  a known  path 
length,  and  then  separated  by  an  electric  field.  Neutral  products  of 
reactions  between  the  two  beams  pass  undeflected  through  the  field  and 
are  subsequently  detected.  When  two  fast  particle  beams  move  parallel 
to  each  other,  not  only  can  their  average  relative  velocity  be  made 
arbitrarily  small,  but  the  velocity  distributions  that  relate  to  the 
inherent  energy  spreads  in  the  beams  axe  dramatically  reduced  as  the 
average  relative  velocity  approaches  zero,  even  though  the  laboratory 
energy  spread  is  unchanged.  Thus  a laboratory  energy  spread  of  1 eV 
in  each  of  two  parallel,  merged  beams  of  equal  mass  and  each  traveling 
at  2.5  keV  (lab)  is  seen  as  a 5 x 10  “*  eV  spread  in  the  center  of  mass 
energies.  The  same  "deamplification"  permits  accurate  control  over 
the  relative  energies  of  the  beams.  Thus,  the  energy  dependence  of 
the  reaction  cross  sections  as  well  as  its  absolute  magnitude  can  be 
established  over  a wide  range  of  energies  (e.g.,  0.15  -*  200  eV), 
providing  more  insight  into  the  details  of  the  reactions  than  is 
possible  by  gas-phase  methods. 

In  practice,  the  low  energy  limit  to  the  merged-beam  cross  section 
measurements  results  from  lack  of  complete  beam  collimation  and 
alignment  and  was  established  to  be  about  0.15  eV.  A practical  limit 
to  the  collimation  results  from  the  need  for  sufficient  ion  beam 
currents  to  attain  useful  signal  levels.  The  method  is  highly 
specific  and  can  provide  great  detail  about  the  reaction  cross  sections. 
It  is,  however,  a complicated  one  and  requires  careful  attention  to 
the  spatial  quality  of  the  bea^s  and  possible  sources  of  interference. 
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Nevertheless,  it  has  been  highly  successful  for  atomic  and  simple 
molecular  ions,  and  has  provided  the  first  real  opportunity  to  test 
theoretical  treatments  of  ion-ion  neutralization.  The  most  serious 
problem  with  the  technique  derives  from  the  fact  that  beams  of 
molecular  ions  (as  comrared  to  atoms)  are  almost  always  at  least 
vibrationally  excitea. 

Little  doubt  exists  as  to  the  correctness  of  the  merged-beams 
data  for  atomic  ions.  The  best-studied  case  is  that  of  H+  + H*  -*  H + H. 
The  low  energy,  merged-beams  data  of  Moseley,  Aberth,  and  Peterson 
(1970)  overlapped  well  the  high-energy  inclined-beam  data  of  Rundel, 
Aitken,  and  Harrison  (1969),  and  detailed  structure  can  be  seen  in 
the  cross  section  (Figure  3).  The  high-energy  data  have  been  improved 
by  Gaily  and  Harrison  (1970),  and  the  intermediate-energy  data  have 
been  improved  by  Peart,  Grey,  and  Dolder  (1976)  who  used  an  inclined 
beam  method.  No  satisfactory,  sophisticated  theoretical  calculation 
is  yet  available  for  H + H . The  problem  is  complicated  by  the 
degeneracies  and  the  need  to  include  rotational  coupling  between 
initial  and  final  states  when  dealing  with  light  molecules. 

The  collected  merged  beams  results  for  atomic  and  molecular 
ions  have  been  presented  by  Moseley,  Olson  and  Peterson  (1975), 
typically  for  the  energy  range  0.15-200  eV.  Reaction  rate  coefficients 
for  use  in  atmospheric  modeling  were  calculated  from  those  data,  using 
a polynomial  fit  (based  on  the  Landau-Zener  theory)  to  extrapolate  to 
lower  energies  and  then  averaging  over  a Maxwellian  velocity 
distribution.  Those  rate  coefficients  are  tabulated  in  Table  1.  A 
direct  comparison  of  the  merged-beams  and  gas-phase  results  will  be 
given  below. 

C.  Flowing  Afterglow  Technique 

At  the  University  of  Birmingham  (UK),  a gas  phase  technique  has 


7 


CROSS  SECTION 


RELATIVE  SPEED  — cm/sec 


Figure  3 


Experimental  data  for  H + H -*  H + H cross 
sections.  Solid  circles:  Moseley,  A berth 
and  Peterson  (1970);  Triangles:  Rundel, 
Altken  and  Harrison  (1969),  and  Gaily  and 
Harrison  (1970);  Continuous  curve  traces  the 
data  of  Peart,  Grey  and  Dolder  (1976). 


Table  I.  ION-ION  MUTUAL  NEUTRALIZATION  RATE  COEFFICIENTS  AT  300  K 
DEDUCED  FROM  MERGED  BEAM  MEASUREMENTS  AT  SRI. 


. -8  3. 

a(10  cm  /sec) 


H + H 

+ - 
N + 0 

+ - 
0+0 

Na+  + O" 

He+  + h" 

He+  + D- 

+ - 
H,  + D 

2+ 

NO  + 0 
+ 

0+0 

2+ 

\ + 0 

N,  + 0, 

2+  2- 
NO  + 0„ 

+ 2- 

°2+  + °2 

N„  + NO  ” 

2+  2- 

NO  + NO, 

+ - 

NO  + OH  • HO 

+ - 2 

0„  + NO 

2+  2- 

NO  + NO, 

+ 3- 

O.  + N0_ 


39  + 21 
12+3  (a) 
20  + 10 
21  + 10 
65  + 15  (b) 
52  + 12  (b) 
47  + 15 
49  + 20 
10+4 
20  + 10 
16+5 
58  + 10 
42  + 13 
13+5 
51  + 15 
10+5  (a) 
41  + 13 
81  + 23 
13+4 


Revised  or  unpublished  data  iron  Bennett  et  al  (1974) 
b + — 

The  data  for  He  + D are  limited  to  energies  2 60  eV  and  do  not 
permit  direct  extrapolation  to  thermal  energies.  However,  the  data 
agree  well  with  theoretical  calculations  that  yield  a 300  K rate 
coefficient  of  52  x 10*  cm3/sec  for  He+  + D and  65  x 10"®  cm3/ sec 


been  developed  for  measuring  ion-ion  mutual  neutralization  rate 
coefficients,  the  Birmingham  method  is  well-suited  for  atmospheric 
related  studies  because  ion-ion  reaction  rate  coefficients  are 
measured  directly,  and  at  atmospheric  temperatures  (measurements  have 
been  made  for  180-530  K).  It  should  also  be  possible  to  measure  three- 
body  reaction  rate  coefficients. 

The  Birmingham  group  uses  a flowing  afterglow  apparatus,  with  a 
movable  Langmuir  probe  to  determine  ion  densities,  and  a mass  spec- 
trometer downstream  to  identify  the  ion  species.  The  afterglow  tube  is 
pyrex,  7.6  cm  in  diameter,  about  84  cm  long.  The  flowing  buffer  gas 
is  usually  helium  at  about  0.5  torr.  Small  concentrations  of  reactant 
gases  may  be  added  either  upstream  or  downstream  of  a microwave 
discharge  cavity  to  create  desired  ionic  species.  Positive  ions  are 
created  either  in  the  microwave  discharge  or  by  Penning  ionization 
in  the  afterglow.  Negative  ions  are  created  by  electron  attachment. 

The  primary  ions  may  be  (and  often  are)  converted  to  different  species 
through  ion-molecule  reactions  taking  place  in  the  flow  tube.  If  the 
ion  densities  are  great  enough,  ion-ion  mutual  neutralization  dominates 
over  ambipolar  diffusion  in  the  loss  of  ionization.  Initial  ion 
densities  in  the  ion-ion  plasma  must  be  typically  3 x 10  cm  to 
avoid  diffusion  effects.  Using  the  Langmuir  probe  to  determine  both 
the  positive  and  negative  ion  densities  as  a function  of  distance 
from  the  ion  production  region,  the  ion-ion  mutual  neutralization  rate 
coefficient  may  be  determined.  The  distance  scale  in  the  flow  tube  is 
converted  to  a time  scale  through  a direct  measurement  of  the  plasma 
flow  velocity,  using  the  Langpiuir  probe  to  follow  a pulsed  disturbance 
down  the  length  of  the  flow  tube. 

Care  must  be  taken  in  the  determination  of  the  ion  identities 
because  the  interpretation  of  the  Lanpnulr  probe  current-voltage 
characteristics  requires  knowledge  of  the  ion  masses,  which  may  be 
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changing  along  the  flow  tube  due  to  ion-molecule  reactions  (Smith  and 
Church,  1976).  In  some  cases,  the  concentrations  of  the  minority 
gases  (which  form  the  ions  in  a neutral  buffer  gas)  may  be  adjusted  so 
that  a single  negative  ion  and  a single  positive  ion  completely 
dominate  the  mass  spectrum.  In  other  cases,  when  one  is  attempting 
to  study  reactions  between  clustered  ions,  a mixture  of  ion  types  is 
likely  to  result  and  some  estimate  of  the  effective  mass  made.  Since 
the  ion-ion  mutual  neutralization  rate  coefficient  may  be  determined 
from  the  decay  of  either  the  positive  or  negative  ion  densities,  the 
investigator  has  some  built-in  duplication  that  aids  in  understanding 
the  ion  chemistry  in  complex  cases. 

Aside  from  the  care  needed  in  assessing  the  ion  chemistry  in  the 
flow  tube,  the  flowing  afterglow  method  appears  to  be  limited  to  the 
study  of  negative  ions  that  have  a fast  attachment  rate  coefficient. 

If  the  electron  attachment  rate  in  the  flow  tube  is  not  large,  an  ion- 
ion  plasma  cannot  be  produced  in  sufficient  density  that  mutual 
neutralization  dominates  over  diffusion.  Often,  species  that  have  a 
large  attachment  rate  also  have  a high  electron  affinity  (EA)  that  tends  to 
lower  the  mutual  neutralization  rate  (Olson,  1972)  though  at  least  one 
known  exception,  SF^  (which  has  a low  EA  but  a high  attachment  rate), 
has  been  studied  in  the  flowing  afterglow  (Church  and  Smith,  1977). 

However,  the  success  of  the  Birmingham  measurements  outweighs 
those  limitations.  Reaction  rate  coefficients  have  been  determined 
for  a large  number  of  ion  species,  both  for  "simple”  and  clustered 
ions,  at  temperatures  relevant  to  atmospheric  studies.  In  Table  2 
we  list  the  data  that  the  Birmingham  group  has  produced.  Where  more 
than  one  ion  species  were  present,  an  average  ion-ion  mutual 
neutralization  rate  coefficient  was  determined. 

Generally,  the  300  K rate  coefficients  determined  with  the  flowing 


Table  II 


ION-ION  MUTUAL  NEUTRALIZATION  RATE  COEFFICIENTS  a MEASURED 
AT  300  K USING  THE  BIRMINGHAM  FLOWING  AFTERGLOW  APPARATUS. 


Ions 


Table  II  Continued 
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a(  10  cm  /sec) 


NH,  (0.16) 

4+ 

NO  (0.16) 


NH,  *NH  (0.23) 

+ 4 3 

NH4  *(NH3)2  (0.26) 
NH4+-NH3  (0.27)1 


N03  (0.19) 


NO  (0.19)  r +-{N0o  *HN0  (0.38) 


5.8  + 0.9 


NH,  • (NH  ) (0.40) 
4+  3 3 

h3o  *(h2o)3  + Cl 
h3o+.(h2o)3  + no3 


NH,  *(NH  ) (0.33)V  + JN0,  -HNO,(0.33) 


N03  *(HN03)2(0.43) 


N03  (0.30) 


N03  .(HN03)2(0.37) 


(and  minority  species) 


H30  + N°3  *i®°3  (an^  minority  species) 

no+-(no2)2  + no3".(hno3)3 


6.1  + 1.0 


4.8  + 0.6 

5.5  + 1.0 
5.7  + 1.0 

3.5  (C) 


Note:  When  more  than  one  positive  or  negative  ion  species  was  present, 
the  fractional  concentrations  are  given  in  parentheses. 

Measurements  at  other  temperatures:  a - 8.3  + 1.0  at  220  K,  and 

+1#3  _g  3 " 

a * 5.6  * at  430  K,  in  units  of  10  cm  /sec. 

k -O  .o 

Measurements  at  other  temperatures:  a - 8.4  + 1.0  at  185  K,  6.4  + 0.8 
at  420  K,  and  5.8  +0.6  at  530  K,  in  units  of  10  ^ cm3/sec. 

CReaction  studied  at  182  K.  If  a ~ T \ then  a • 4.5  x 10~®  cm^/sec  at 
182  K can  be  used  to  derive  a rate  coefficient  at  300  K as  given  above, 


afterglow  apparatus  are  lower  than  those  deduced  from  the  merged-beams 

data,  although  there  is  little  direct  overlap  (two  cases)  between  the 

ion  species  studied.  Both  sets  of  results  tend  to  be  within  a factor 

of  2 of  the  absorbing-sphere  theoretical  values,  with  the  flowing 

afterglow  results  lower  and  the  merged  beams  results  higher.  Mutual 

neutralization  rate  coefficients  for  atomic  ions  show  a large  variation, 

which  is  understandable  in  terms  of  the  varied  locations  of  potential 

curve  crossings,  the  number  of  crossings,  and  the  ion  masses.  The 

flowing  afterglow  data  for  molecular  ions  seem  to  imply  that  the 

presence  of  many  electronic  states  and  of  internal  degrees  of  freedom 

tends  to  narrow  the  range  of  values  of  the  rate  coefficients,  although 

the  effects  of  mass  and  electron  affinity  are  certainly  present.  The 

data  obtained  on  the  Birmingham  apparatus  thus  far  yield  rate  coef- 

•8  3 

ficients  in  the  range  4-10  x 10  cm  /sec  and  Indicate  that  a value 
-8  3 

of  6 x 10  cm  /sec  at  300  K is  a good  estimate  of  the  mutual 
neutralization  rate  coefficient  for  the  complex  molecular  ions  if  one 
wishes  to  choose  a single  value,  without  regard  to  the  precise  species 
involved. 

We  have  stated  earlier  that  the  flowing  afterglow  rate  coef- 
ficients (300  K)  tend  to  be  smaller  than  those  deduced  from  the  merged- 

beams  data.  Only  two  reactions  have  been  studied  in  common  between 

■f  • + • 

those  apparatuses,  NO  + NO^  and  NO  + NO^  . Unfortunately,  the  merged- 
beams  results  for  those  two  reactions  were  much  larger  than  the  body 
of  results  for  all  other  reactions  studied  with  that  apparatus,  and 
the  belief  is  that  those  beam  ions  were  partially  excited.  However, 
it  is  still  interesting  to  see  a comparison  between  the  flowing 
afterglow  data  and  the  merged-beams  data  for  a particular  case. 

No"*"  + NO ^ . For  this  comparison,  the  flowing  afterglow  thermal  rate 
coefficients  were  converted  to  monoenergetic  rates  (cross  section 


times  relative  velocity)  and  the  results  are  shown  in  Figure  4.  Two 
points  are  made:  the  merged-beams  results  are  almost  10  times  greater 
than  the  flowing  afterglow  results  (and  the  two  cannot  be  reconciled 
by  a different  extrapolation  procedure),  and  the  energy  dependence  of 
the  flowing  afterglow  data  agrees  with  theoretical  expectations  at 
low  energies,  Qv  ~ E whereas  the  merged-beams  data  show  a slower 
decrease  with  energy. 


MERGED-BEANS  DATA, 
EXTRAPOLATED  USING 

£ + B + Cv  + Dv2 


FLOWING  AFTERGLOW  DATA, 
EXTRAPOLATED  USING 


0.1  1 
RELATIVE  ENERGY  — eV 


Figure  4 A comparison  of  merged-beams  data  and  flowing  afterglow 
data  for  the  ion-ion  mutual  neutralization  reaction 
NO  + NO-  -•  neutral  products 


As  a test  for  excitation  in  the  merged-beam  ions,  one  series  of 
experiments  was  made  to  determine  if  the  NO**"  ions  were  in  the 
metastable  a3£+  state.  Charge  transfer  attenuation  of  N0+  beams 
passing  through  Ar  gas  were  made  in  the  manner  of  Turner,  Rutherford 


and  Compton  (1968),  and  it  was  concluded  that  less  than  10%  of  the 

ions  were  excited.  Similar  tests  were  made  to  see  if  NO  ions 

2 

formed  in  different  gases  were  vibrationally  excited,  but  the  tests 
were  inconclusive.  Very  likely,  both  the  N02  and  NO^  ions  were  at 
least  vibrationally  excited,  and  the  possibility  also  exists  that 
some  of  the  NO^  ions  were  in  the  isomeric  "peroxy"  state.  Because 
the  ground  electronic  states  of  the  negative  ions  are  generally 
"offset"  geometrically  from  the  ground  neutral  states,  vibrational 
excitation  could  enhance  the  Franck-Condon  factors  for  transitions 
between  these  states,  and  the  neutralization  rates  for  vibrationally 
excited  ions  could  be  substantially  larger  than  those  for  "cold"  ions. 
Similarly,  the  transition  probabilities  could  be  larger  for  the 
peroxy  N0^  than  for  the  normal  form.  Thus,  the  results  from  the 
Birmingham  apparatus  are  more  trustworthy. 


; 


Ill  THEORY 


The  basic  ideas  in  the  theoretical  description  of  ion-ion  mutual 
neutralization  have  been  reviewed  by  Moseley  et  al.  (1975).  With  the 
assumption  that  electron  transfer  takes  place  at  avoided  crossings  of  the 
ionic  and  covalent  potential  curves  for  the  interacting  systems,  the 
Landau-Zener  theory  may  be  used  to  calculate  the  mutual  neutralization 
cross  section  provided  that  (a)  the  coupling  matrix  element  can  be  calcu- 
lated for  each  crossing,  (b)  the  excited  states  of  the  neutral  system  are 
known  in  order  to  determine  the  asymptotic  covalent  potential  curves,  and 
(c)  not  too  many  states  are  involved  so  that  the  problem  may  be  tractable. 
An  appropriate  formula  for  the  matrix  elements  connecting  the  initial  and 
final  states  has  been  deduced  by  Olson,  Smith,  and  Bauer  (1971)  by  corre- 
lating a large  number  of  ab  initio  calculations  and  experimental  data  on 
electron  transfer.  In  regard  to  the  number  of  covalent  states  Involved, 
there  are  two  extremes: 

(a)  If  only  one  or  two  covalent  states  are  possible  as 
product  channels,  a close-coupling  calculation  may 
be  carried  out  instead  of  using  the  Landau-Zener 
theory.  The  close-coupling  calculation  takes  account 
of  interactions  over  an  extended  region  of  internu- 
clear  separations  near  an  avoided  crossing  rather  than 
just  at  the  point  of  closest  approach  between  two 
interacting  states.  Olson  (1977)  has  carried  out 
close-coupling  calculations  for  Na+  + Cl"  and  1C*-  + C1", 
where  there  is  nominally  only  one  avoided  crossing, 
but  spin-orbit  splitting  induces  a second  crossing. 

For  KC1 , it  is  noteworthy  that  Olson  found  Landau- 
Zener  cross  sections  to  be  several  orders  of  magnitude 
too  small  (compared  with  the  more  reliable  close- 
coupling results)  because  the  avoided  crossings  take 
place  at  large  internuclear  distance  (~  20  £)  where 
the  Landau-Zener  description  of  the  interaction  is 
inadequate. 
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If  there  is  a large  number  10)  of  covalent  states 
crossing  the  Ionic  potential,  as  is  common  for 
molecular  systems,  Olson  (1972)  has  developed  an 
"absorbing  sphere"  model  that  avoids  calculating 
transition  probabilities  for  each  curve  crossing. 

The  model  seems  to  be  quite  successful,  and  the 
results  are  dependent  only  on  the  electron  detach- 
ment energy  of  the  negative  ion  and  the  reduced  mass 
of  the  system  (Moseley  et  al.,  1975). 


Emphasis  in  recent  years  has  been  on  mutual  neutralization  of  large 
ionic  clusters,  motivated  by  the  realization  that  such  clusters  dominate 
in  the  nighttime  quiescent  atmosphere  at  altitudes  below  about  85  km  when 
negative  ions  are  most  important.  With  molecular  ions,  one  interesting 
problem  is  to  determine  how  excitation  of  internal  vibrational  and  rota- 
tional degrees  of  freedom  during  the  collisions  affects  the  mutual 
neutralization  cross  sections.  If  the  excitation  takes  place  at  long 
ranges,  the  cross  section  may  be  increased,  because  absorption  of  the 
available  kinetic  energy  (kT)  can  trap  the  ions  in  orbiting  states  and 
increase  the  probability  of  electron  transfer.  The  clustering  of  neutral 

molecules  to  an  ion  core,  e.g.,  H 0+,(H  0)  , will  affect  the  mutual 

3 2 n 

neutralization  cross  sections  both  through  altered  curve  crossings  (both 
in  number  and  crossing  distance)  and  by  increasing  the  number  of  internal 
degrees  of  freedom.  In  addition,  the  increased  mass  of  clustered  ions 
will  tend  to  lower  the  reaction  rates  because  of  the  decrease  in  thermal 
velocities.  The  SRI  group  has  considered  the  fact  that  the  bonding  of 
enough  water  molecules  to  an  ion  could  result  in  the  lowering  of  the  ionic 
curve  below  all  covalent  ones,  so  that  mutual  neutralization  could  occur 
only  by  coalescence  of  the  positive  and  negative  ions  into  an  ionic 
complex,  with  Internal  (vibrational)  excitation  absorbing  the  binding 
energy  (Bennett  et  al.,  1974). 

A detailed  calculation  of  absolute  cross  sections  for  the  mutual 
neutralization  of  large  clustered  ions  is  extremely  difficult  because  of 
the  large  number  of  electronic  states  and  external  degrees  of  freedom. 


A more  reasonable  goal  is  to  determine  the  dependence  of  the  cross 
sections  on  the  number  of  added  molecules,  their  dipole  moments,  and 
their  binding  energies. 

We  will  describe  here  the  historical  developments  in  our  research 
to  this  problem.  In  an  earlier  report  we  gave  estimates  of  neutralization 
rates  for  hydrated  ions  (Bennett  et  al.,  1974).  We  have  examined  the 
ion-ion  neutralization  theory  for  hydrated  (or  other  cluster)  ions.  This 
extends  our  earlier  results,  which  applied  to  moderately  large  clusters 
(with  average  hydration  number  n > 6) , to  smaller  clusters  (1  < n < 6) . 

Olson's  absorbing  sphere  version  of  the  curve-crossing  theory  of 
neutralization  is  highly  effective  for  many  small  molecular  ions,  but  it 
is  probably  inapplicable  for  the  hydrated  (cluster)  ions.  The  absorbing 
sphere  model  does  not  consider  that  the  initial  ion  and  product  neutral 
minimum  energy  states  probably  have  substantially  different  configurations 
and  intemuclear  separations,  with  correspondingly  small  Franck-Condon 
factors  between  the  initial  and  final  states  of  the  electron  capture 
transitions.  Such  an  effect  would  reduce  the  transition  probabilities. 
Likewise,  the  model  does  not  account  for  the  "tidal"  orbital  trapping 
that  can  result  from  the  excitation  of  internal  vibrational  and  rotational 
modes  during  the  collision.  That  effect  will  increase  the  transition 
probabilities.  To  consider  these  effects,  we  first  took  a radically 
different  approach  to  obtain  upper  and  lower  limits  to  neutralization 
rates,  independent  of  electronic  transitions.  Part  of  it  follows  the 
results  of  our  earlier  considerations  (Bennett  et  al.,  1974). 

For  cluster  ions  of  any  size,  close  encounters  within  a distance  r 

c 

of  typical  molecular  dimensions  will  lead  to  a hard  impact,  with  almost 
certain  internal  excitation  or  fragmentation,  permanent  trapping,  and 
inevitable  neutralization.  The  term  r depends  on  n but  not  on  the 
relative  energy  E^.  The  cross  section  depends  on  and  fi  as 


o — 


ne2r  (n) 
c 

E 

r 


(3) 
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This  hard-impact  mechanism  leads  to  a lower -bound  rate  coefficient  for 
hydrated  ions  of  magnitude 


for  any  n £ 1. 

Collisions  outside  the  close-encounter  range  r but  inside  some 
larger  tidal  trapping  range  rfc  may  also  lead  to  neutralization,  but  not 
necessarily.  The  tidal  trapping  range  r^  depends  on  as  well  as  on  n, 
and  the  trapping  probability  may  be  much  less  than  unity,  especially  for 
small  n (1  S 5 S 5);  the  trapping  cross  section  is  then  less  than  some 
upper  bound  a , 


After  trapping,  neutralization  follows  with  almost  unit  probability  if 

there  are  enough  internal  degrees  of  freedom,  i.e.,  if  fl  is  moderate  to 

large  3);  it  may  be  less  probable,  however,  for  n near  unity. 

If  we  write  k for  the  rate  coefficient  for  neutralization  by  the 
n,t 

tidal  trapping  mechanism,  the  total  neutralization  rate  coefficient  is 
the  sum: 


The  total  k is  less  than  an  upper  bound  k related  to  o . k (T.n)  4 
n u u n 

^(1,5).  On  the  basis  of  a simple  electrostatic  model,  we  previously 
estimated  k to  be 


W*'*'*‘e  expected  to  be  a good  approximation  to  the  close- 

coupling encounter  rate  k , the  difference  (k  -k  ) is  not  necessarily  a 


good  estimate  of  the  tidal  neutralization  rate  k , and  may  be  much 

n»  t 

larger.  That  is  because  k is  an  estimate  based  on  a static  distortion 

u 

model  that  neglects  dynamic  effects.  Especially  when  the  number  of 

internal  degrees  of  freedom  is  small,  as  when  fi  = 1 or  2,  these  dynamic 

effects  may  severely  limit  the  probability  of  tidal  trapping  or  of 

subsequent  neutralization.  Consequently,  we  expect  k to  be  small 

n , t 

when  n is  small  (perhaps  even  smaller  than  k^) , and  then  to  rise  gradually 

with  increasing  n to  some  limit  approaching  (k  -k  ) . Correspondingly, 

u X 

the  total  neutralization  rate  will  rise  from  near  k.  at  low  n to  near  k 

J l u 

at  large  n. 

This  general  picture  seems  to  be  in  reasonable  agreement  with  the 

experimental  information,  especially  the  data  on  small  cluster  ions  from 

the  Birmingham  group,  which  generally  show  room- temperature  rates  near 
-8  3 

6 x 10  cm  /sec  (Smith,  Church,  and  Miller,  1978). 

A.  P.  Hickman  (1978)  contributed  the  next  step  in  the  theoretical 
treatment,  developing  a model  for  the  quantitative  calculation  of  the 
0 + N0+  and  0 + NO^H^O  reaction  rates.  The  model  considers  separately 

both  the  electron  transfer  and  internal  excitation  processes.  Only  the 
lowest  frequency  rotational  and  vibrational  modes  were  considered,  but 
the  effective  increase  in  the  ionization  energy  of  the  clustered  NO+  ion 
was  included.  The  effects  of  orbiting  on  the  electron  transfer  probabil- 
ities were  not  examined,  but  otherwise,  the  effects  of  internal  excitation 
were  found  to  be  minimal  in  such  simple  clusters;  the  dominant  effects  of 
clustering  were  caused  by  changes  in  the  mass  and  in  the  electron 
affinity.  An  approximate  scaling  formula  was  obtained  that  fits  a wide 
range  of  experimental  data  involving  simple  ions  and  small  clusters. 

This  scaling  rule  relates  the  rate  constant  for  recombination  (cr)  to  the 
electron  affinity  of  the  electron  donor  (EA)  and  the  reduced  collision 
mass  (m) : 

a'1  = C(T/300  K)0*5  m°’5  (EA)°‘4 


If  the  units  of  <*,  m,  and  EA  are  cm  /sec,  atomic  units  and  eV, 

4 

respectively,  then  C ■ 4.38  x 10  . This  formula  fits  a wide  variety  of 
experimental  data  with  an  accuracy  of  about  + 30%,  as  shown  in  Figure  S 
A preprint  of  a paper  describing  this  work  is  included  as  an  Appendix. 


ION  - ION  NEUTRALIZATION  RATE 
Experiment 
$ Simple  Ions 
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Figure  5.  Comparison  between  recent  theoretical  results  and  experi 
mental  data  from  several  laboratories . The  complex  pote; 
approach  does  not  invoke  adjustable  parameters.  Units  a 
follows:  a:  10®  cm3/ sec;  m:  atomic  units;  EA:  electron 


IV  ATMOSPHERIC  IMPLICATIONS 


The  process  of  mutual  neutralization  of  ions  in  the  atmosphere  is  a 
complicated  function  of  altitude.  Three  altitude-dependent  variables 
must  be  considered;  (1)  the  type  of  positive  and  negative  ions  present; 

(2)  the  temperature;  and  (3)  the  neutral  gas  density.  In  the  stratosphere, 
below  about  30  km  altitude,  the  ions  are  all  in  the  form  of  large  clusters . 
Since  the  stratosphere  contains  few  free  electrons,  mutual  neutralization 
is  the  dominant  process  for  the  removal  of  ions.  Furthermore,  the  gas 
density  is  high  enough  that  three-body  neutralization  is  effective.  At 
higher  altitudes,  the  ion  species  become  somewhat  simpler,  and  the  «as 
density  is  such  that  binary  mutual  neutralization  dominates.  A good 
summary  of  average  mutual  neutralization  rate  coefficients  in  the  atmosphere 
appears  in  a chart  by  Smith  and  Church  (1977),  which  is  reproduced  as 
Figure  6.  The  figure  shows  likely  positive  and  negative  ion  types  and  the 
appropriate  ionic  recombination  coefficients  deduced  from  flowing  after- 
glow and  other  experiments. 


The  average  rate  coefficients  are  not  based  on  measurements  for  all  known 
atmospheric  ion  species,  convoluted  with  an  atmospheric  model.  Rather, 
they  are  based  on  the  fact  that  the  binary  rate  coefficients  measured  with 
the  Birmingham  flowing  afterglow  apparatus  do  not  vary  much  with  ion 
identity  at  a fixed  temperature,  and  that  the  temperature  dependence  of 
the  measured  rate  coefficients  seems  to  fit  the  T ^ dependence  predicted 
by  theory.  The  three-body  estimates  are  perhaps  cruder.  Smith  et  al. 
(1978)  conclude  that  the  temperature  dependence  of  the  average  mutual 
neutralization  rate  coefficient  is  probably  a more  important  consideration 
in  atmospheric  modeling  than  is  precise  knowledge  of  the  ion  species. 

In  the  disturbed  atmosphere,  the  ionospheric  layers  are  lowered  in 
altitude.  Depending  on  the  type  of  disturbance  (solar  proton  shower, 
nuclear  blast,  etc.)  there  may  be  a temperature  change  associated  with 
the  altitude  adjustment,  and  there  will  also  be  a trend  toward  simpler 
ion  species  at  a given  altitude,  as  the  ionization  source  strength 
increases.  If  the  amount  of  ionization  versus  altitude  can  be  specified 
for  a certain  disturbance,  then  Figure  6 can  be  adjusted  to  give  approxi- 
mate mutual  neutralization  rate  coefficients  for  the  disturbed  atmosphere. 
The  only  radical  change  that  might  be  expected  would  be  at  high  altitudes 
(>  70  km),  where  a significant  concentration  of  atomic  and  unclustered 
molecular  ions  could  lead  to  a wide  range  of  mutual  neutralization  rate 
coefficients,  depending  on  the  ion  types  present.  In  general,  the  simpler 
atmospheric  ions  should  lead  to  higher  neutralization  rates  (e.g.: 

3 x 10  ^ cm^/sec  for  + O and  0+  + 0 ). 


V FUTURE  WORK 


A great  deal  has  been  accomplished  thus  far  in  understanding  ion-ion 
mutual  neutralization.  The  basic  electron  transfer  mechanism  is  clear. 


For  atomic  ions,  good  calculations  have  been  shown  to  be  feasible  when 
the  computation  time  is  warranted  (Olson,  1977) . Approximate  theory  has 
been  worked  out  for  molecular  ions  (Olson,  1972),  and  the  first  steps 
have  been  taken  to  treat  cluster  ions  quantitatively  (Hickman,  1978). 
Theoretical  research  is  in  progress  at  SRI  to  determine  how  those  ideas 
must  be  modified  for  larger  ions,  the  clusters  that  occur  in  the  atmos- 
phere. Experimentally,  merged-beams  and  inclined-beams  apparatus  have 
been  used  to  measure  mutual  neutralization  reaction  cross  sections  over  a 


wide  energy  range,  and  basic  curve-crossing  theoretical  models  can  account 
for  the  general  characteristics.  Flowing-afterglow  experiments  have 
provided  reaction  rate  coefficients  for  simple  and  clustered  atmospheric 
molecular  ions,  at  atmospheric  temperatures. 

The  thrust  of  future  work  should  be  to  reach  a better  understanding 
of  the  neutralization  of  clustered  ions  in  the  atmosphere.  Theoretical 
research  can  be  expected  to  determine  cross-section  dependencies  on 
molecular  properties  such  as  the  degree  of  clustering,  bond  strengths, 
and  dipole  moments.  Many  more  results  from  flowing-afterglow  experiments 
on  ion  cluster  neutralization  will  be  necessary  to  verify  the  trends 
predicted  by  theory.  Furthermore,  such  data  would  test  the  preliminary 
indication  from  flow  tube  results  that  the  neutralization  rate  coefficients 
do  not  vary  much  (a  factor  of  2)  between  different  clusters.  The  flowing 
afterglow  method  can  also  be  used  to  study  three-body  mutual  neutralization 
and  presumably  such  work  will  be  done. 

Technology  exists  today  that,  in  principle,  would  allow  merged-beams 
experiments  to  be  carried  out  with  clustered  ions,  using  high-pressure. 


r 


nozzle-expansion  beam  sources  and  cryogenic  vacuum  pumping.  However, 

such  experiments  are  very  difficult  and  the  rewards  do  not  currently  i ] 

11 

seem  worth  the  effort  and  expense,  in  view  of  declining  budgets. 

Finally,  we  should  point  out  that  knowledge  of  the  cross  section  is 
a solution  to  only  half  of  the  ion-ion  mutual  neutralization  problem. 

• ! I 

Very  little  is  known  about  the  final  states  of  the  neutral  products  (see 
Moseley  et  al.,  1975).  However,  workers  at  the  University  of  Birmingham 
have  recently  used  a monochromator  with  their  flowing  afterglow  apparatus 
to  observe  radiation  resulting  from  ion-ion  mutual  neutralization 
reactions. 


A« 
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APPENDICES 


T.  M.  Miller  of  SRI  spent  the  period  16  May-30  September  1977  working 
on  ion-ion  mutual  neutralization  experiments  at  the  University  of 
Birmingham,  using  a flowing-afterglow  apparatus  with  Langmuir  probe 
diagnostics.  A paper  was  published  recently  in  the  Journal  of  Chemical 
Physics  on  the  flowing-afterglow  experiments  at  Birmingham  during  that 
period  of  collaboration,  and  its  abstract  is  attached  here  as  Appendix  A. 

In  the  University  of  Birmingham  Ionic  Physics  Laboratory,  the  flowing 
afterglow  apparatus  shares  some  crucial  equipment  with  a Selected  Ion 
Flow  Tube  (SIFT)  apparatus.  During  Dr.  Miller's  stay  in  Birmingham,  he 
took  part  in  SIFT  experiments  in  which  ion-molecule  reactions  were  studied 
in  regard  to  the  stratosphere.  The  abstract  of  the  paper  describing  that 
work  appears  here  as  Appendix  B. 

Appendices  C through  F are  abstracts  of  the  theoretical  papers  already 
published  under  this  contract,  and  Appendix  G is  a preprint  of  a theoretical 
paper  to  be  published  in  1979  in  the  Journal  of  Chemical  Physics. 


Appendix  A 

J.  Chem.  Phys.  68,  1224  (1978) 


Mutual  neutralization  of  simple  and  clustered  positive  and 
negative  ionsa) 

David  Smith,  Michael  J.  Church,  and  Thomas  M.  MiHeH1* 

Department  of  Space  Research,  University  of  Birmingham,  Birmingham  BIS  2TT.  England 
(Received  14  September  1977) 

Measurements  are  reported  of  the  rate  coefficients,  a,  for  several  ion-ion  mutual  neutralization  reactions 
principally  involving  NH,*  ions  and  their  ammonia  clusters  NH^<NHj)uj  with  several  different  negative 
ions.  The  data  were  obtained  utilizing  an  ion-ion  flowing  afterglow  plasma  combined  with  Langmuir 
probe  diagnostics.  Moat  of  the  measurements  were  obtained  at  300  K although  the  NH<  +CT  reaction 
has  also  been  studied  at  220  and  430  K.  Both  the  absolute  magnitude  of  a and  its  temperature  variation 
are  shown  to  be  in  acceptable  agreement  with  theoretical  predictions.  The  a for  the  cluster  ion  reactions 
are  very  similar  to  those  for  the  simpler  ions,  even  when  both  positive  and  negative  ions  are  large  dusters, 
although  an  increase  in  the  mean  reaction  cross  section  for  the  duster  ion  reactions  is  discernible.  All  the 
a measured  to  date  for  both  simple  and  duster  ion  reactions,  albeit  for  spades  of  high  electron  affinity 
and  over  the  limited  temperature  range  of  110  to  330  K,  are' within  the  range  (4-IO)x  10~*  cm1  a*1,  more 
■mHh«i  variations  occurring  with  temperature  than  with  ionic  mass. 
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J.  Chem.  Phys.  69,  308  (1978) 


A laboratory  study  of  the  reactions  of  N+,  Na+,  N,+, 
N4+,  0+,  02+,  and  NO+  ions  with  several  molecules  at 
300  K 


D.  Smith,  N.  G.  Adams,  and  T.  M.  Miller*)  b> 

Department  of  Space  Research.  University  of  Birmingham.  Birmingham  B15  2TT,  England 
(Received  17  January  1978) 

A study  has  been  made  of  the  rate  coefficients  and  product  ion  distributions  for  the  reactions  at  300 K of 
the  ions  N\  N,\  N,\  N.\  O*.  O,*.  and  NO*  with  CH,NH„  NH„  H,S,  CH,OH,  H,CO.  COS,  O* 
H]0,  CH„  COj,  CO,  H,,  and  N,  molecules  listed  in  tncreasini  order  of  their  kmizatioa  energies.  These 
measurements  are  intended  as  a contribution  to  stratospheric  chemistry.  In  the  binary  reactions  of  the 
ions  of  large  recombination  energy  with  molecules  of  low  ionization  energy,  multiple  ion  products 
generally  result  and  the  rate  coefficients  are  dose  to  gas  kinetic.  Conversely,  the  low  recombination 
energy  ions  NO*  and  Of  generally  undergo  ternary  association  reactions  with  the  large  ionisation  energy 
molecules.  The  reactions  of  NJ  and  N,*  are  very  similar,  the  most  common  mechanism  apparently  being 
direct  charge  transfer  usually  followed  by  fragmentation,  the  nitrogen-nitrogen  bonds  in  the  reacting  ions 
remaining  intact.  The  N*  and  N,*  reactions  differ  foom  the  N,  and  N4*  reactions  in  that  they  show  a 
greater  propensity  to  form  N-X  bonds,  X«0.  C S,  H.  etc.  The  O*  and  Of  reactions  generally 
proceed  via  direct  charge  transfer  where  energetically  poeefolt 


Appendix  C 


f 


■ 

Phys.  Rev.  Lett.  35,  1073  (1975) 


Semiclassical  Perturbation  Scattering  by  a Rigid  Dipole* 

F.  T.  Smith,  D.  L.  Hues!  is,  and  D.  Mukhorjee 
Stanford  Research  Institute,  Menlo  I'art;,  California  91025 

and 

W.  H.  Miller 

Department  of  Chemistry,  University  of  California,  Hcrkeley,  California  9-1720 

(Received  HI  duly  1075) 

A uniform  semi  classical  S matrix  has  liccn  developed  for  collisions  of  charged  parti- 
cles with  rotating  rigid  dipoles,  with  use  of  first-order  perturbation  theory.  The  re- 
sulting expression  is  analytical,  depending  on  tabulated  functions,  and  trivial  to  calcu- 
late; it  allows  evaluations  of  quantum  transitions  in  classically  forbidden  regions,  and 
of  quantum  interference  effects. 
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Phys.  Rev.  A 17 , 939  (1978) 


Semiclassical  perturbation  theory  of  electron-molecule  collisions 

W,  H.  Miller 

Department  of  Chemistry.  University  of  California.  Berkeley.  California  94720 

F.  T.  Smith 

SRI  International.  Menlo  Park,  California  94025 
(Received  25  July  1977) 

The  theory  of  the  classical  or  semiclassical  S matrix  is  combined  with  the  use  of  perturbation  dynamics  to 
derive  an  approximately  unitary  expression  for  the  scattering  matrix  for  a very  general  class  of  potential 
interactions.  The  S matrix  takes  the  form  of  a sum  over  products  of  licsset  functions  whose  orders  are 
related  to  the  changes  in  quantum  numbers  occurring  in  the  transition,  and  whose  arguments  depend  on  the 
dynamical  variables  of  the  problem,  including  the  unperturbed  quantum  numbers.  In  the  general  case,  these 
arguments  can  be  expressed  as  simple  integrals  over  the  unperturbed  trajectory,  and  for  well-behaved 
potentials  they  can  be  explicitly  evaluated  in  terms  of  the  modified  Bessel  functions  JC„  and  K,.  The 
connection  between  the  semiclassical  perturbation  scattering  theory  and  other  approximations,  such  as  the 
Born  and  eikonal  approximations,  is  demonstrated  The  general  theory  is  illustrated  by  applications  to 
electron-  (or  ion-)  polar-molecule  scattering,  including  quadruple  as  well  as  dipole  interactions  and 
including  coupling  to  vibrjlions  in  both  harmonic  and  enharmonic  approximations  I he  more  complicated 
interactions  involve  lengthier  products  of  Bessel  functions  in  the  siim-and-product  representation,  but  these 
are  easily  and  systematically  evaluated,  and  they  reduce  smoothly  to  the  appropriate  simpler  expressions 
when  the  coupling  coefficients  of  the  higher-order  terms  become  small  More  complicated  potentials, 
including  interactions  between  polyatomic  molecules,  can  be  handled  by  a simple  systematic  extension  of  the 
same  principles.  For  electron-molecule  scattering,  these  expressions  can  be  used  in  their  present  form  since 
the  sums  are  dominated  by  Bessel  functions  of  comparatively  low  order  which  can  be  evaluated  directly; 
extensions  to  molecule-molecule  scattering  and  ion-molecule  scattering  are  equally  valid  formally,  but  their 
practical  application  will  offen  require  the  use  of  asymptotic  approximations  to  the  Bessel  functions. 
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Phys.  Rev.  A L7,  954  (1978) 


Semielassieal  perturbation  theory  of  electron-polar-molecule  collisions 
Total  excitation  and  scattering  cross  sections 


D.  Mukherjec*  and  F.  T.  Smith 

Molecular  Physics  Center.  SRI  International,  Menlo  Park.  California  9 4025 
(Received  25  July  1977) 

Semielassieal  perturbation  scattering  theory  is  applied  to  electron -polar-molecule  scattering  Cross  sections 
for  elastic  and  rotationally  inelastic  transitions  arc  obtained  for  Ay  0,  I,  2.  3,  4,  in  the  case  of  a pure 
charge-dipole  interaction  potential.  Results  are  presented  in  terms  of  functions  of  dimensionless  parameters 
involving  the  moment  of  inertia  / and  dipole  moment  of  the  target,  and  the  mass,  charge,  and  energy  of 
the  projectile.  A range  of  these  panimeters  sufficient  to  describe  most  situations  of  practical  interest  is 
explored.  Significant  oscillations  arc  found  in  the  cross  sections  for  |Ay  | --  0 and  2 as  a function  of  dipole 
strength  /J  - mpe/li*.  For  |Ay'|  = I,  the  Dorn  approximation  is  shown  to  be  appropriate  for  /3<  I,  and 
invalid  for  /J  ~ 1 or  >1  Scaling  rules  are  deduced  which  should  aid  in  the  correlation  and  extrapolation  of 
quantal  calculations  on  specific  systems.  With  the  guidance  of  (he  Born  approximation,  scaling  principles  are 
also  suggested  for  differential  scattering. 
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Phys.  Rev.  A 17,  968  (1978) 


Momentum  transfer  in  electron-polar-molecule  collisions:  Results  of 
semiclassieal  perturbation  scattering  theory 

A.  P.  Hickman  and  F.  T.  Smith 
Molecular  Physics  Center,  SRI  International.  Menlo  Park,  California  94025 
(Received  25  July  1977;  revised  manuscript  received  l.f  February  1978) 

Semiclassieal  perturbation  scattering  theory  is  applied  to  electron  polar-molecule  collisions,  assuming  a 
pure  charge-dipole  interaction  potential.  Cross  sections  for  momentum  transfer  are  given  in  terms  of 
dimensionless  parameters  involving  the  moment  of  inertia  and  dipole  moment  of  the  target,  and  the  mass, 
charge,  and  energy  of  the  projectile.  A range  of  these  parameters  ts  explored  corresponding  to  collisions  at 
energies  above  0.1  cV  for  most  molecules.  For  large  values  of  the  di|to|e  moment,  the  momentum  transfer  is 
much  smaller  than  the  prediction  of  the  first-order  Born  approximation,  and  involves  significant  contributions 
both  from  elastic  scattering  and  from  transitions  with  large  values  of  JA^J.  The  behavior  of  the  cross 
section  as  a function  of  (he  rotational  quantum  number  j is  found  to  obey  a simple  scaling  rule 


approximate  scaling  formula  that  fits  a wide  range  of  experimental  data  to 
an  accuracy  of  about  + 30%. 


MP  78-114 
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I Introduction 


This  paper  concerns  the  theoretical  study  of  the  Mutual  neutralisation 


reaction  between  positive  and  negative  ions.  This  general  class  of  reactions 


where  X and  Y are  charged  species  that  suiy  range  in  else  and  complexity 


from  simple  ions  such  as  H , 0 , 0, 


or  NO  to  hydrated  or  "cluster1 


The  hydration  or  cluster  number 


ions  such  as  HO 


n stay  range  from  1 to  6 or  even  larger.  Typically  one  of  the  neutral 


products  of  the  reaction  is  in  an  excited  state,  denoted  in  Eq.  (1)  by 


X . This  class  of  reactions  is  responsible  for  the  final  removal  of  ions 


at  altitudes  below  80  km  in  the  ionosphere,  and  knowledge  of  the  rate 


constants  is  of  considerable  practical  isq>ortance 


The  purpose  of  the  present  work  is  to  develop  a model  that  can  be 


used  to  estismte  rate  constants  for  a wide  variety  of  reactions  of  the 


form  (1).  TWo  distinct  mechanisms  will  be  incorporated  into  the  model 


The  first  is  neutralisation  by  electron  transfer,  which  has  been 


previously  studied  by  Olson.  Our  treatsmnt  of  electron  transfer  is 


similar  to  Olson's.  However,  the  present  theory  is  more  general  because 


it  provides  a unified  framework  within  which  both  electron  transfer  and 


internal  excitation  can  be  treated.  Internal  excitation  has  been  proposed' 


as  a mechanism  leading  to  a capture  reaction  and  formation  of  neutral  XY 


The  interaction  of  the  projectile  with  the  target  causes  a transfer  of 


energy  AE  from  the  initial  translational  mode  to  internal  vibrational 
and  rotational  modes.  If  AE  is  greater  than  the  initial  relative  kinetic 
energy,  the  collision  partners  will  be  trapped  in  a large  ellipsoidal 

orbit.  If  energy  transfer  back  to  translational  motion  is  unlikely,  as 
one  would  expect  when  the  target  is  a large  cluster,  then  over  a period 
of  several  orbits,  each  successive  close  encounter  will  lead  to  greater 
internal  excitation  and  a more  tightly  bound  orbit.  The  ultimate  result 
will  be  a larger,  neutral  cluster. 

In  Section  II  we  develop  in  detail  the  theoretical  model,  which  is 

based  on  treating  the  loss  of  ions  to  the  neutralization  channels  by 

adding  an  imaginary  part  to  the  potential.  The  connection  with  Olson's 

absorbing  sphere  model  is  discussed.  Numerical  results  are  presented 

and  discussed  in  Section  III.  Our  calculations  Indicate  that  electron 

transfer  is  the  dominant  neutralization  mechanism  for  collisions  involving 

small  ions  and  clustered  species  with  n < 3.  The  calculations  also  yield 

scaling  rules  for  such  collisions.  We  illustrate  the  effect  of  single 

hydration  by  comparing  in  detail  collisions  of  0 + N0+  and  0~  + NO+*H  0, 

2 

and  discuss  qualitatively  the  application  of  the  model  to  larger  clusters. 
Section  IV  contains  concluding  remarks. 


II  Theory 

A.  Complex  Potential  Model 

In  this  section  we  develop  a complex  potential  model  to  treat  ion-ion 
mutual  neutralisation.  We  first  consider  the  case  of  the  electron  transfer 
mechanism,  and  then  show  how  Internal  excitations  can  alao  be  Included. 
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A schematic  diagram  of  the  relevant  potential  curves  is  shown  in  Figure  1. 

The  Idea  Is  that  charge  exchange  can  occur  via  a curve  crossing  from  the 

initial  ionic  potential  to  any  one  of  a number  of  closely  spaced  states  of 
* 

the  form  X + Y.  It  Is  assumed  that  the  final  statea  are  sufficiently  closely 

spaced  that  at  any  separation  R of  X+  and  Y* , some  final  state  curve  will  cross 

the  initial  state.  At  this  R,  the  transition  occurs  predominantly  to  this 

particular  final  state.  The  matrix  element  connecting  the  initial  and  final 

statea  will  be  denoted  by  Hfi(R).  An  approximate  formula  for  Hfi(R)  has  been 

3 

deduced  by  Olson,  Smith,  and  Bauer  by  correlating  a large  number  of  ab  initio 

calculations  and  experimental  results. 

Our  model  consists  of  treating  the  closely  spaced  final  states  as  a 

continuum,  and  defining  an  R-dependent  transition  probability  or  width  T. 

The  function  T(R)  gives  the  rate  of  transitions  out  of  the  initial  ionic 

channel  at  each  R.  The  dynamics  of  the  collision  are  treated  by  adding  an 

imaginary  component  - “ T(R)  to  the  potential  function  V(R)  of  the  ionic 

channel.  (In  the  present  case,  V(R)  • - 1/R).  This  complex  potential  model 

4-6 

has  been  used  in  various  other  applications  involving  an  initial  state 
embedded  in  a continuum  of  states,  and  its  application  to  an  initial  state 
crossing  a series  of  discrete  but  closely-spaced  states  has  been  discussed  by 
Miller  and  Morgner . ^ The  general  conditions  necessary  for  the  validity  of  the 
s»del  are  the  following:  1)  The  final  statas  must  be  sufficiently  closely 
spaced  that  transitions  may  be  considered  to  occur  continuously  along  a 
trajectory,  rather  than  at  Isolated  curve  crossings;  2)  the  probability  of 
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transitions  from  the  neutral  channels  back  to  the  Ionic  channel  must  be 


negligible.  Ue  expect  these  conditions  to  be  valid  in  the  present  case. 

We  now  formulate  the  model  more  precisely.  We  define  (in  the  standard 

way  ) 


r(R)  - 2ttp|h£1(r)|2 


(2) 


Pis  the  density  of  final  states.  For  a manifold  of  final  states  whose  energy 
levels  are  given  (in  a.u.)  by  a Rydberg  formula 


E 

n 


pis  easily  shown  from  the  work  of  Miller  and  Morgner7  to  be 


(3) 


P - [ 2(EA  + R-1) J (*) 

where  EA  is  the  electron  affinity  of  Y.  We  adopt  Eq.  (4)  as  the  definition 
of  p even  though  the  electronic  energy  levels  of  a complicated  molecule  may 
deviate  somewhat  from  Eq.  (3). 

The  matrix  element  H££(R)  is  defined  by  Olson,  Smith,  and  Bauer3  and 
Olson*  as  follows: 

Hfi(R)  ■ 1.044  If*  R*  exp  (-0.857  R*)  (5) 

where 

R*  • <0  + (6) 

and 


The  quantity  1^  ia  the  detachment  energy  (in  atomic  units)  of  the  electron 
being  transferred.  For  ground  state  negative  ions,  I is  just  the  electron 


affinity  EA.  Formula  (5)  is  estimated  by  Olson  to  be  valid  to  about  a 


The  probability  of  electron  transfer  on  a trajectory  whose  orbital 
angular  momentum  is  i is  given  by* 


where 


In  this  formula  E is  the  collision  energy,  and  m is  the  reduced  mass.  For 


trajectories  in  the  Coulomb  potential,  the  classical  turning  point  (distance 
of  closest  approach)  R is  given  by 


R - (1  / m)(l  + (1+2E1  /m),J 


The  total  cross  section  for  neutralisation  is  given  by 


rz  TX&: ar.' 


where  k is  related  to  the  collision  energy  E by  E ■ A /(2m)  and  to  the 
velocity  v by  v - ftk/m.  Ihe  rate  constant  a is  obtained  by  taking  the 
appropriate  thermal  average  over  a Boltzmann  distribution  of  velocities 


f(v,T): 


where 


<*(T)  - J*  f(v,T)  a(v)  v3  dv 
0 


For  thermal  energies,  essentially  all  of  the  velocity  dependence  of  the 

2 

cross  section  Eq.  (12)  is  contained  in  the  1/k  prefactor.  This  is  because 
bhe  Coulomb  field  accelerates  particles  to  velocities  many  times  greater 
than  their  initial  thermal  velocities.  To  a good  approximation,  therefore, 

in  E<1*  ^ may  ke  considered  independent  of  energy,  and  R in  Eq.  (11) 

X 

may  be  taken  to  be 


R * l /2m 
X 


Then  if  v ■ (8kTAnm)  ' is  the  average  velocity  at  temperature  T,  it  is  easily 


shown  that 


er(T)  - “ v <j(v) 


We  now  consider  Che  extension  of  Che  shove  model  Co  include  neuCrslisstion 


vis  internal  excitation  of  the  target.  It  is  necessary  to  include  terms  in 
the  potential  corresponding  to  the  internal  vibration  and  rotation  of  the 
target.  We  consider  here  the  case  that  the  target  has  a dipole  moment  and 
a single  vibrational  coordinate.  We  will  treat  the  rotational  and  vibrational 

excitation  using  the  formulas  of  Semiclasslcal  Perturbation  Scattering  (SPS) 

v 9-11 
theory. 

In  the  absence  of  the  electron  transfer  mechanism,  the  SPS  formulas 

10 

presented  by  Miller  and  Smith  may  be  used  to  calculate  the  transition 
probabilities  from  the  initial  state  to  Che  final  state  Jj^i^v^. 

These  are  obtained  from  the  T matrix  elements,  which  are  written  as  simple 
analytic  functions  of  the  "average"  quantum  numbers  J * ^(j.  + jf),  etc., 
and  the  changes  AJ  ■ (J^  - j^),  etc.: 


T - T(J ; J I v;  Aj  LI  Av)  (17) 

Kukherjee  and  Smith11  showed  that  the  total  cross  section  for  a particular 
(j^)  -•  (jfv^)  transition  is  given  by 


a(Jivi  - Vi> 


-T— 2 Z 12  (2J  + 1)|T(J;J  A v;  AJ  Li  Av)  |2  (18) 

It  (2Jt  + l)  I Ai  J 

where  the  sum  is  over  values  of  J,  Li,  and  J consistent  with  angular  moamntum 
conservation.  (Explicit  limits  are  given  by  Mukherjee  and  Smith.11) 


Formula  (18)  may  be  written  in  the  form 


°(jiVi  ~ jfVf)  " ^2  2 (M  + DQj  ^ivi  “*  J£v£>  (19) 

k 7 

where 

QI  ^iVi  "*  ^fvf*  “ (2j  + 1) (21+  1)  2 2(2J  + 1>IT<J*  Tv;  Aj  A1Av)|2  (20) 

1 J 

The  similarity  of  Eq.  (19)  to  Eq.  (11)  is  evident.  If  one  identified  7 with 
the  angular  momentum  l in  the  spherically  symmetric  case,  one  may  interpret 
Q£  in  Eq.  (19)  as  the  probability  for  a particular  degeneracy  averaged) 
transition  on  a trajectory  whose  average  orbital  angular  momentum  is  7. 

The  preceding  formulas  for  Q7  and  o(j  v - j v,)  are  based  on  the 
assumption  that  electron  transfer  cannot  occur.  If  we  relax  this  assumption, 
then  a single  collision  may  lead  either  to  electron  transfer,  or  to  Internal 
excitation  without  electron  transfer.  The  cross  section  for  the  former 
process  is  given  as  before  by  Eq. (12).  We  define  the  cross  section  for  the 
latter  process  by 


1 

1 


®(Vi  - jfV  “ *2  2(2J+1)  vVi  - jfV  (2l) 

k i 

where 

Q/Vi  - J£vf>  - (1  - p,)  - V£>  (22) 
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In  ocher  words,  Che  probability  of  inCernal  excitation  on  a given 


Crajeccory  depnds  on  Che  probability  Chat  electron  transfer  does  not  occur. 

This  assumption  is  completely  analogous  to  the  assumption  frequently  made 

in  treating  spherically  symmetric  complex  potentials,  that  the  real  part  of 

the  phase  shift  does  not  change  as  the  imaginary  part  of  the  potential  is 
12 

"turned  on."  Such  a perturbation  assumption  is  consistent  with  our  use 
of  a perturbation  approach  to  the  internal  excitation. 

Formulas  (12)  and  (21)  may  be  used  to  calculate  cross  sections  for 
electron  transfer  and  internal  excitation  of  an  ionic  polar  target.  In 
order  to  calculate  the  contribution  to  the  total  cross  section  for  neutralization 
it  is  necessary  to  make  the  further  assumption  that  an  internal  excitation 
energy  AE  leads  invariably  to  recombination  when  AE  > E (E  is  the  initial 
translational  energy).  This  is  an  excellent  approximation  as  long  as  there 
are  sufficient  Internal  modes  to  make  an  energy  transfer  back  to  the 


i 


El 

■'t'  * 


translational  mode  unlikely. 

B.  Relation  to  Absorbing  Sphere  Model 

The  present  complex  potential  approach  may  be  viewed  as  a generalization 

of  the  absorbing  sphere  model^  (ASM).  The  ASM  assumes  that  if  the  collision 

partners  approach  closer  than  a certain  critical  distance  , then  charge 

transfer  via  curve  crossing  occurs  with  unit  probability.  R depends  on  the 

c 

effective  mess  m and  electron  affinity  EA,  and  is  the  distance  at  which  the 
matrix  element  for  curve  crossing  first  reaches  s particular  threshold  value. 

The  cross  section  for  neutrellsatlon  is  given  by 


-x1 ' ’ ; •••  . ■ > 
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I 


^ — r 


iKT 

m 

■r . 


ttR 

o - + ttR 

L C 


(23) 


At  thermal  energies  this  Is  well  approximated  by 


nR 


a e= 


E (2*> 

The  complex  potential  model  assigns  a neutralization  rate  (transition 
probability  per  unit  time)  to  each  value  of  the  coordinate  R.  The  formulas 
(8) -(12)  correspond  to  integrating  this  transition  probability  along  the 
classical  trajectories  determined  by  the  potential  V(R)  “ - 1/R.  Formula 
(11)  relating  the  distance  of  closest  approach  and  the  angular  momentum  of 
a Coulomb  trajectory  enables  us  to  compare  the  results  typically  obtained 
for  the  complex  potential  model  for  [Eq.  (8)]  with  the  corresponding 
prediction  of  the  absorbing  sphere  model.  The  comparison  is  shown 


schematically  in  Fig.  2.  The  ASM  gives  a step  function  whose  discontinuity 


is  at  l ■ (2mR  ) , whereas  the  complex  potential  model  tends  to  smmoth  out 


the  discontinuity. 

The  use  of  an  R-dependent  transition  probability  makes  the  present  model 
somewhat  more  flexible  than  the  ASM.  This  flexibility  may  be  especially 
important  when  one  or  both  of  the  ionic  species  are  hydrated.  The  qualitative 


discussion  in  Section  III-B  shows  that  the  energetics  of  cluster  formation 


imply  that  f(R)  should  be  zero  for  values  of  R less  than  some  R*,  where  R* 


will  depend  on  the  size  of  the  cluster.  For  sufficiently  large  clusters,  R 
may  be  comparable  to  R^.  In  this  case,  the  probability  of  neutralization  by 
the  curve  crossing  mechanism  would  be  drastically  reduced.  For  such  a large 
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cluster,  one  might  expect  the  internal  excitation  mechanism  to  be  much  more 
important.  The  complex  potential  model,  with  a r related  to  cluster  size, 
provides  a consistent  framework  for  examining  the  collision  process  as  a 
function  of  cluster  size. 


S': 


III  Results  and  Discussion 

A.  Predicted  Scaling  Behavior 

We  have  found  that  the  results  of  the  complex  potential  model  for 
a wide  variety  of  systems  can  be  summarized  by  a simple  scaling  rule.  Here 
we  discuss  how  this  rule  was  obtained. 

It  was  first  determined  that  the  electron  transfer  mechanism  was 
much  more  important  than  internal  excitation  for  simple  ions  and  small 
clusters.  When  only  the  former  mechanism  is  included,  the  numerical 
calculation  of  the  rate  constant  « for  any  particular  system  depends  only 
on  m and  EA.  We  then  sought  to  express  the  results  of  all  the  numerical 
calculations  by  a formula  of  the  form 


Or(T)  « (T/300  K)‘°‘5  m"°*5  <EA)X 


(25) 


•0.5 


The  T ‘ behavior  arises  whenever  the  cross  section  varies  as  v‘2.  This 
behavior  was  obtained  by  Olson  in  the  ASM,  and  was  verified  in  the  present 
case  by  direct  numerical  calculation.  The  m"0*5  behavior  was  also  observed 
in  the  calculations.  It  is  easily  understood  by  returning  to  the  ASM, 
which  predicts 


_"*c 

or  ~ v o(V)  - v — — 
fi 
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(26) 


•»«**  if 


t-.  T • . 


_2  — 0.5 

Since  E ~ %mv  and  v * (8kT/nm)  * , 


0.5  0.5 

m T 


(27) 


However,  our  examination  of  the  absorbing  sphere  model  has  shown  that  R 

c 

depends  most  strongly  on  EA,  and  only  logarithmically  on  m.  We  therefore 

expect  that  the  dependence  of  q on  m and  EA  might  be  approximated  by  a 

power  law  as  in  Eq.  (25).  If  this  formula  is  reasonable,  then  a log-log 
0.5 

plot  of  a*J>  («t  T • 300  K)  vs  EA  will  give  a straight  line  of  slope  \. 
Fig.  3 shows  that  such  a plot  is  nearly  linear.  The  slope  is  determined  by 
a least  squares  fit  to  be  X - -0.4.  We  therefore  obtain  the  following 
formula 


a"1  - C(T/300  K)0,5  m°'5  (EA)°*4  (28) 

3 

If  the  units  of  or,  m,  and  EA  are  cm  /sec,  au.,  and  eV,  respectively, 
then  C - 4.38  x 104. 

Figure  4 and  Table  II  show  that  Eq.  (28)  provides  a simple  correlation  of 
13-17 

wide  range  of  data.  Considering  that  the  basic  formula  for  Hfi(R)  was  only 

expected  to  be  accurate  to  a factor  of  three,  and  that  the  data  range  from 

^ jm 

collisions  of  H and  H to  H^O  *(^2^3  an<*  ^3  * the  approximate  scaling 
rule  must  be  considered  remarkably  successful. 


B.  Qualitative  Discussion  of  the  Effects  of  Clustering 

The  energetics  of  cluster  formation  have  been  previously  discussed,2 
and  in  some  cases  quantitative  calculations  have  been  carried  out.^*^ 

This  section  discusses  the  application  of  the  complex  potential  model  to 
larger  clusters. 
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In  general,  one  expects  that  the  addition  of  successive  waters  of 

hydration  will  lower  the  energy  of  the  ion  pair  state  relative  to  the 

2 

corresponding  neutrals.  The  situation  is  illustrated  by  the  schematic 
diagram  in  Fig.  5.  The  energy  of  X+,H20  is  lower  than  that  of  X+  and  an 
infinitely  separated  H20  because  of  the  charge-dipole  and  polarization 

if 

interactions.  In  contrast,  one  expects  the  energy  of  X ‘I^O,  at  the 

equilibrium  separation  of  X+*H20,  to  be  slightly  higher  than  the  energy  of 

* 

X . 

* 

The  curve  X + Y in  Fig.  5 denotes  the  lowest  state  to  which  charge 

transfer  can  occur.  Because  of  the  stabilizing  effect  of  clustering,  discussed 

in  the  preceding  paragraph,  this  asymptote  moves  closer  and  closer  to  the  ionic  state 

as  the  cluster  size  increases.  In  the  context  of  the  complex  potential 

model,  we  may  interpret  this  to  mean  that  T(R)  should  be  nonzero  only  for 
* 

R > R . An  estimate  previously  used  is  that  each  water  of  hydration  moves 
the  two  asymptotes  1 eV  closer  together.  If  this  is  the  case,  then  R*  is 
given  (in  a.u.)  by 

1/R*  * W - n/27.2  (29) 

n is  the  hydration  number,  and  W is  the  energy  difference  (in  a.u.)  between 
, - + * 

the  asymptotes  X + Y and  X + Y , and  may  be  crudely  estimated  by 

V fet  IP(Y)  - EA(X)  (30) 

that  is,  the  difference  between  the  ionisation  potential  of  Y 
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and  Che  electron  affinity  of  X.  For  many  species  of  atmospheric  interest, 

W is  in  the  range  6-8  eV  (0.2-0. 3 a.u.). 

Trial  calculations  were  performed  in  which  T(R)  was  set  to  zero 
★ 

for  R < R >10  aQ.  The  change  in  the  rate  constant  was  negligible  (less 
than  17.),  indicating  that  electron  transfer  occurs  mainly  at  large  R.  It 
is  likely  that  cluster  formation  would  cause  a reduction  in  the  matrix 
element  for  electron  transfer.  We  found  that  reducing  r(R)  by  a factor  of 
four  caused  only  about  a 257.  decrease  in  the  rate  constant.  Combined  with 
Eq.  (29)  and  our  estimate  of  W,  these  results  suggest  that  for  small 
clusters  (n  ;£  3),  hydration  does  not  significantly  inhibit  electron  transfer. 


SSi 


C.  The  Internal  Excitation  Mechanism:  Calculations  for  0 + NO 


and  0 + NO  -HjO 


The  system  0 + NO  is  particularly  well  suited  to  the  application 
of  the  complex  potential  model  because  the  real  part  of  the  potential,  both 


the  isotropic  and  anisotropic  parts,  can  be  modelled  asymptotically  using  the 


20  + 

dipole  moment  function  of  NO  . Cross  sections  for  rotational  and 


vibrational  excitation  have  been  calculated  both  with  and  without  the 
imaginary  part  of  the  potential,  using  Semiclasslcal  Perturbation  Scattering 
(SPS)  theory.  The  application  of  this  theory  in  the  case  where  the  dipole 


moment  is  a linear  function  of  the  vibrational  coordinate  has  been  discussed 


10 


by  Miller  and  Smith.  Furthermore,  we  have  contructed  a simple  model  for 


the  geometry  of  NO  -H^O,  and  have  examined  the  cross  sections  for  excitation 


of  the  lowest  frequency  rotational  and  vibrational  modes.  The  results  allow 
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us  to  examine  the  competition  between  the  electron  transfer  and  the  internal 


excitation  mechanisms  for  simple  ions,  and  also  to  obtain  information  about 


the  changes  that  occur  with  the  first  hydration 


Figure  6 shows  the  assumed  geometry  of  the  NO  *11,0  complex.  For 


simplicity,  we  assume  that  the  center  of  charge  is  at  the  midpoint  of  the 


and  that  the  interaction  of  the  oxygen  of  the  water  molecule  with  the 


N and  the  other  0 can  be  described  by  a Born -Mayer  potential  of  the  form 


where  A ■ 74.445  a.u.  and  B ■ 2.006  a.u.  This  potential  is  the  geometric 


mean  of  the  Born-Mayer  potentials  for  neutral  0-N  and  0-0  interactions 


The  lnternuclear  separation  of  NO  is  taken  to  be  the  unperturbed  value 


2a  ■ 2.00  a • With  these  assumptions , the  equilibrium  geometry  will  have 


0 ■ tt/2,  and  R can  be  determined  from  the  minimum  value  of  the  Interaction 


where  u ■ 1.85  D ■ 0.728  a.u.  is  the  dipole  moment  of  H_0.  For  these 


values  of  the  parameters,  the  eqiullbrium  value  of  R is  4.9  a.u.,  and  the 


The  above  model  allows  us  to  calculate  the  dipole  moment  (relative  to 


t of  inertia,  and  rotational  energy  level  spacing 


of  NO  *10,  treating  it  as  a rigid  linear  molecule.  These  values  are 


auanarised  in  Table  1.  The  audel  also  allows  us  to  estimate  the  frequencies 
of  the  lowest  vibrational  or  bending  modes,  and  to  obtain  the  dependence  of 


the  dipole  moment  function  on  the  corresponding  internal  coordinates.  We 
consider  both  vibration  in  the  R coordinate  and  bending  in  the  9 coordinate. 
Estimated  frequencies  for  these  modes,  assuming  they  are  uncoupled,  are  also 
given  in  Table  I. 

The  results  of  our  calculations  may  be  sunmarized  as  follows.  For 
0 + N0+,  the  electron  transfer  mechanism  dominates  the  scattering.  Any 
trajectory  that  would  lead  to  significant  rotational  or  vibrational  excitation 
in  the  absence  of  the  width  r,  has  a high  probability  (essentially  unity)  of 
electron  transfer  when  T is  turned  on.  Also,  for  those  trajectories  where 
the  electron  transfer  probability  is  small  (i.e.,  large  l or  large  classical 
turning  point),  vibrational  excitation  is  much  less  important  than  rotational 
excitation  and  may  be  neglected.  We  have  also  found  that  electron  transfer 
is  the  most  important  mechanism  for  0 + NC^H^O,  but  that  Internal  excitation 

may  increase  the  rate  constant  by  about  107.  at  300  K.  The  reason  for  this 
is  that  the  cluster  has  a larger  dipole  moment  than  the  bare  NO*,  and  hence 
rotational  excitation  can  occur  for  more  distant  trajectories  than  for  NO* 
alone.  The  calculations  showed  that  the  cross  sections  for  vibrational 
excitation  were  considerably  smaller  than  those  for  rotational  excitation 
for  all  trajectories  not  leading  to  electron  transfer.  In  other  words,  for 
the  purposes  of  the  complex  potential  model,  the  singly  hydrated  cluster 
behaves  as  if  it  were  a rigid  rotor. 

The  final  rate  constants  obtained  were  as  follows.  For  0 + N0+, 

-18  3 

or  (300  R)  ■ 13.9  x 10  cm  /sec,  essentially  all  from  electron  transfer. 


5' 


For  0 + N0+ *H20,  the  rate  constant  from  electron  transfer  is  a (300  K)  - 
-8  3 

13.1  x 10  cm  /sec;  the  inclusion  of  internal  excitation  changes  this  to 

14.1  x 10  ® cm^/sec. 

The  extension  of  the  present  model,  with  its  detailed  assumptions 
about  geometry,  to  larger  clusters  is  probably  not  justified  without 
rigorous  structure  calculations.  A few  general  consents  may  be  pertinent, 
however.  The  addition  of  more  water  molecules  will  probably  decrease,  not 
increase,  the  total  dipole  moment,  since  some  cancellation  could  be  expected 
to  occur.  This  suggests  that  large  rotational  excitation  due  to  distant 
trajectories  may  be  important  only  for  single  (or  un-)hydrated  ions.  The 
large  clusters  will  have  a greater  number  of  internal  modes  that  can  be 
excited  by  direct  impact.  Previous  discussion  has  indicated  that  electron 
transfer  may  not  occur  at  close  distances  for  these  larger  clusters.  It 
may  be  that  at  small  distances  neutralization  occurs  instead  with  high 
probability  by  the  internal  excitation  mechanism. 

IV  Conclusions 

He  have  presented  a theory  of  ion-ion  mutual  recoefeinatlon  based  on 
modelling  the  loss  of  particles  from  the  initial  channel  to  a manifold  of 
final  channel  by  adding  an  imaginary  part  to  tha  Coulomb  potential.  We 
have  examined  the  relative  Importance  of  competing  mechanisms,  and  have 
found  that  the  electron  transfer  process  accounts  for  at  least  SOX  of  the 
rate  constant  in  the  cases  considered.  The  results  of  the  theory  have 
been  parameterized,  and  provide  a useful  and  reasonably  accurate  correlation 
of  a wide  range  of  experimental  data. 


This  correlation  formula,  Eq.  (28),  predicts  a small  decrease  in  the 
neutralization  rate  as  the  size  of  the  positive  ion  cluster  increases.  This 
effect  is  due  entirely  to  the  change  in  the  effective  mass  of  the  relative 
motion.  However,  our  more  detailed  study  of  collisions  of  hydrated  N0+  with 
0 showed  a net  increase  in  the  rate  as  the  hydration  number  increased  from 
zero  to  one,  because  of  the  increased  probability  of  internal  excitation. 

It  is  felt  that  these  results  cannot  be  interpreted  to  indicate  any  general 
trend  of  neutralization  rates  as  a function  of  cluster  size.  They  appear 
only  to  suggest  that  the  rates  for  small  clusters  are  not  substantially 
different  from  the  rates  for  the  corresponding  bare  ions. 
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Table  I.  Values  of  the  molecular 


in  the  calculations.  Note  that  the  vibrational  coordinate  of  N0+  is 


the  N-0  distance;  that  of  NO  -HjO  is  the  coordinate  R shown  in  Pig.  6 
with  NO  at  its  equilibrium  value. 


Rotational  constant  (cm  ) 


Moment  of  inertia  (gm  cm  ) 


Equilibrium  value  of  vibrational 
coordinate,  R (a  ) 


Equilibrium  value  of  dipole 
relative  to  c.m. , u,  • u 


Dipole  derivative  (R  ) (a.u.) 


Vibrational  Frequency  (a.u.)  w 


Bending  Frequency  ( 


n 

Table  II.  Details  of  the  experimental  data 
plotted  in  Figure  4.  Units  are  as  in  Fig.  4. 


Ions 

m (EA) 

or 

Reference 

H+  + h" 

27 

39  ± 21 

13 

h2+  + d' 

38 

47  ± 15 

17 

N-+  + 0 " 

2 2 

118 

16  ± 5 

14 

H+  + o" 

137 

26  ± 8 

13 

o+  + o" 

142 

27  ± 13 

13 

Na+  + o" 

153 

21  ± 10 

15 

N0+  + o" 

162 

49  ± 20 

15 

o2+  + o' 

164 

10  ± 4 

15 

NO+  + N02‘ 

232 

6.4  ± 0.7 

16 

NH,+  + Cl" 

4 

247 

6.7  ± 0.7 

16 

nh4+.(nh3)2  + no2" 

270 

4.9  ± 0.6 

16 

NO+  + N03’ 

278 

5.7  ± 0.6 

16 

nh4+.(nh3)2  + Cl' 

329 

7.9  ± 1.0 

16 

H20+-(H20)3  + Cl' 

350 

4.8  ± 0.6 

16 

h3o+  *(h20)3  + no3' 

356 

5.5  ± 1.0 

16 

CC1F2+  + Cl' 

359 

4.1  A 0.4 

16 

CC13+  + Cl" 

375 

4.5  ± 0.5 

16 
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Figure  Captions 

1.  Schematic  diagram  of  the  relevant  potential  curves  for  ion-ion  mutual 
neutralization.  The  state  X*  + Y is  one  of  a manifold  of  closely  spaced 
states  indicated  by  the  light  lines. 

2.  Schematic  illustration  of  the  neutralisation  probability  as  s 

function  of  orbital  angular  momentum  i or  classical  turning  point 
2 

R - l /2m,  for  the  complex  potential  swdel  (solid  line)  and  the 

X 

absorbing  sphere  model  (dashed  line). 

3.  Demonstration  that  the  numerical  results  of  the  complex  potential 

model  (points)  may  be  approximately  fit  by  a straight  line  on  an 

appropriate  log-log  plot,  leading  to  the  power  law  dependence  of 

Eqs.  (25)  and  (28).  EA  is  given  in  eV,  m in  a.u.,  and  or  in 
-8  3 

10  cm  /sec. 

4.  Comparison  of  the  scaling  law,  Eq.  (28),  and  experimental  data. 

or  is  calculated  at  T ■ 300  K.  EA  is  given  in  eV,  m in  a.u.,  and 
*8  3 

a in  10  cm  /sec.  The  data  is  tabulated  for  reference  in  Table  II. 

5.  Schematic  diagram  showing  how  hydration  tends  to  lower  the  energy  of 
an  ion  pair  state  relative  to  the  corresponding  neutrals.  In  the  text 
it  is  argued  that  for  R < R , T(R)  should  be  sero. 

6.  The  geometry  assumed  for  a simple  model  of  NO* .I^O.  The  diagram 
is  not  to  scale. 
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